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Abstract 

Coordinated multi-point (CoMP) communication is attractive for heterogeneous cellular networks 
(HCNs) for interference reduction. However, previous approaches to CoMP face two major hurdles 
in HCNs. First, they usually ignore the inter-cell overhead messaging delay, although it results in 
an irreducible performance bound. Second, they consider the grid or Wyner model for base station 
locations, which is not appropriate for HCN BS locations which are numerous and haphazard. Even 
for conventional macrocell networks without overlaid small cells, SINR results are not tractable in the 
grid model nor accurate in the Wyner model. To overcome these hurdles, we develop a novel analytical 
framework which includes the impact of overhead delay for CoMP evaluation in HCNs. This framework 
can be used for a class of CoMP schemes without user data sharing. As an example, we apply it to 
downlink CoMP zero-forcing beamforming (ZFBF), and see significant divergence from previous work. 
For example, we show that CoMP ZFBF does not increase throughput when the overhead channel delay 
is larger than 60% of the channel coherence time. We also find that, in most cases, coordinating with 
only one other cell is nearly optimum for downlink CoMP ZFBF. 



I. Introduction 

Improving cellular network capacity is a serious concern of network operators, as new genera- 
tion mobile devices (e.g. smart phones and mobile connected tablets/laptops) keep getting wider 
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adoption and generating crushing data demand Jill. One approach among many viable solutions 
is coordinated multi-point communication, where multiple cells cooperate to improve the key 
quality-of-service metrics including network throughput (see [2] for an overview). As the cellular 
network is now in a major transition from conventional macrocell network to heterogeneous 
cellular network with additional tiers of small cells (e.g. microcells, picocells, femtocells and 
distributed antennas), it is highly desirable to study CoMP techniques in this new paradigm. 

A. CoMP Study Hurdles in Heterogeneous Cellular Networks 

Existing research studies the CoMP concept in conventional macrocell-only networks [121-0. 
However, two important aspects of existing works prevent their direct application to the new 
environment of HCNs. 

The first one is the idealized assumption on overhead messaging among neighboring cells, 
i.e. assuming no inter-cell overhead delay. Not surprisingly, it results in promising predictions 
on CoMP performance such as multi-fold throughput gain [|2l-[l8]|. However, inter-cell overhead 
delay is not trivial in typical cellular network jQl-JTSl, which leads to an irreducible performance 
bound in theory [14] and significant performance degradations in practice [l9ll- [fTT1l . The impor- 
tance of inter-cell overhead delay is further confirmed by latest industrial implementations, where 
the performance degradations are much smaller when inter-cell overhead channel is particularly 
optimized (e.g. directly connecting base stations with gigabyte Ethernet) ifTOl . ifTTI . Clearly, 
inter-cell overhead delay is an important performance limiting factor, and must be modelled and 
quantified in CoMP study ifTSl . However, this is not a trivial task in HCNs where different types 
of base stations (BSs) have very different backhaul capabilities and protocols ifTTI - lfTBl . [fT6ll . 

The second one is the SINR characterization. Previous works use the grid model or the Wyner 
model of base station locations to characterize the end-user's signal and other-cell interference 
in CoMP schemes. Unfortunately, neither model is suitable for HCNs because they assume 
base stations (BSs) are located on regular positions (e.g. locations form a hexagon [HI, ifTTl . a 
line [[TSl or a circle [7J) while small cells in HCNs have unplanned ad hoc locations. Besides, 
SINR characterization under these two models are inaccurate or intractable, or both. The Wyner 
model allows clean-form SINR results, to understand CoMP concept from information-theoretic 
perspective. However, this model is not accurate due to unrealistic assumptions on the wireless 
channel and inter-cell interference [19J. On the other hand, the grid model of conventional 
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macrocell networks is known to be intractable for SINR analysis (H. In HCNs with additional 
tiers of overlaid small cells, the grid model is more complex and tractable SINR results are 
almost hopeless [20], not to mention a grid model for small cells is unlikely to be very realistic. 
In sum, previous models are incapable to capture the new characteristics of BS locations in 
HCNs, and new models are highly desired for accurate and tractable SINR analysis. 

B. Previous Work 

Early theoretical literature completely ignores the impact of overhead messaging in CoMP 
schemes, i.e. they assume that overhead messages have no quantization error and the overhead 
channel is delay-free with infinitely large capacity [|2l-[[6l. Such an ideal assumption is useful 
for the understanding of CoMP fundamentals, but is obviously far from reality in most practical 
cases. As a result, it causes highly over-optimistic predictions on the performance of CoMP 
schemes [l9l-lfTll. 

Practical issues of overhead messaging are considered in a few more recent works. The 
capacity limit of inter-cell overhead channel is considered especially in CoMP joint processing 
where user data is shared among cells ETI . [|22l . The limited feedback model is widely used 
to characterize the quantization inaccuracy in overhead messages ifTTl . [|23l . [[24|. However, the 
impact of overhead delay is either ignored or considered under a very simplified model - that is, 
a fixed delay model [25] . In general, appropriated modeling of overhead delay is still missing, to 
capture the imperfections of overhead channel such as congestion and hardware delays. In our 
previous work [|26l. we proposed various models on overhead channels in HCNs (e.g. backhaul 
and over-the-air overhead channels) and derived the respective delay distributions. These results 
will be used in this paper to quantify the impact of overhead delay in CoMP performance. 

Because the grid model and the Wyner model are obviously not suitable for SINR charac- 
terization in HCNs, several recent works focus on developing new models for BS locations 
[|20l . [|27l . These works model the locations of BSs in HCNs as nodes in one or more spatial 
Poisson Point Processes (PPPs). Base station transceiver parameters (e.g. transmit power and 
path-loss exponent) become the mark of the respective node in the PPP. In this way, the PPP 
model characterizes the BS location randomness as well as the heterogeneity among different 
types of BSs. Previous studies on this PPP model show that, besides providing analytical 
tractability, the new model is at least as accurate as the hexagon-grid model in characterizing the 



4 



SINK distribution [|20||, [|271|-[|29l|. Therefore, some in industry have begun to use it for SINK 
characterization in HCNs [l30ll . 

C. Contributions 

This paper evaluates downlink CoMP in HCNs, using appropriate models of inter-cell overhead 
delay and BS locations. We first develop a new analytical framework for the evaluation of a class 
of CoMP schemes without user data sharing among cells. This framework quantifies the impact 
of inter-cell overhead delay in HCNs by using our previous results on the delay distributions 
under various overhead channel configurations [|26l|. Note that this framework includes previous 
CoMP analysis without overhead delay modeling as a special case. Therefore it can be used to 
explain the performance gaps between previous analytical predictions and real implementations 
of CoMP 

To concretely illustrate the usage of this framework, we apple it to a specific scheme: downlink 
CoMP zero-forcing beamforming (ZFBF). CoMP ZFBF has been studied in macrocell networks 
before |[7l, [fT7|. [[3T| and is attracting industrial implementation efforts ^H^. We derive upper and 
lower bounds on the end-user SIR distribution for CoMP ZFBF in HCNs, using the spatial PPP 
model to characterize other-cell interference from all BSs in the entire plane. These bounds are 
closed-form and show clear dependence on important parameters such as the overhead message 
bit size and the number of coordinated cells. Using this SIR characterization along with the 
CoMP evaluation framework, we quantify the downlink CoMP ZFBF coverage and throughput 
as functions of overhead messaging configurations. Compared with previous work, our results 
provide new design insights for CoMP ZFBF, for example, on the best number of coordinated 
cells and the appropriate configuration of overhead channels. 

II. System Model 

A. Downlink Heterogeneous Cellular Network Modeling 

We consider a downlink heterogeneous cellular network consisting of K different types of 
base stations (e.g. macrocells, microcells, picocells, femtocells and distributed antennas). We 
refer to a specific type of BSs as a tier, and thus call the network a J^-tier HCN. For the A;-th 
tier, the BSs have transmit power P^., number of antennas N^, path loss exponent ak and spatial 
density BSs per unit area. For example, compared with macrocells, femtocells typically have 
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much lower transmit power, fewer antennas and eventually a much higher spatial density as tens 
to hundreds of femtocells will often be deployed in the area of a macrocell [fT6l . 

We consider a typical end-user equipped with a single antenna. We denote its location as the 
origin and the locations of BSs as {Xi^, k = 1,2 . . . , K, i E N}, where is the location of 
the i*^ closest BS to the origin in the /c*'^ tier. We assume all tiers are independently distributed 
on the plane and BSs in the /c*'^ tier are distributed according to homogeneous Poisson Point 
Process (PPP) $fc with intensity A^. 

For the purpose of cell association, the end-user will listen to the downlink pilot signals from 
different BSs, and measure their long-term average powers. With short-term fading averaged 
out, the end-user will associate with the BS from whom it receives the strongest average power 
max jPfclXj fcl""''}. The index of the selected serving BS is 

fc=l,2...,X,ieN I . I J 

ii\k*} = aig max |PJXi J'^n = arg max = 1, PJXi J""*}. (1) 

Therefore, we denote the selected serving BS as BSi fc*. When BSi fc* is capable of coordinating 
with L other cells, it usually would prefer to coordinate with the L strongest interfering cells, 
for example, to maximally cancel inter-cell interference in CoMP zero-forcing beamforming [l2l, 
ifTTl . We denote the set of these L coordinated BSs as Bl- 

B. Overhead Messaging in CoMP Schemes 

In this paper, we investigate the impact of realistic inter-cell overhead signaling on a class of 
CoMP schemes, where an end-user is served by only one BS without user data sharing among 
coordinated cells. To cooperate with the serving cell, a coordinated BS BSj a: G Bl needs to be 
updated with certain key parameters in that cell. The choices of these cooperation-dependent 
parameters are different among various CoMP schemes, with common examples being user 
channel states and user scheduling information. These parameters naturally fluctuate over time 
because of the dynamics in network environment and user mobility. 

Assumption 1: (i.i.d. discrete time model) We assume the cooperation-dependent parameters 
stay constant for a time Ti^k (either deterministic or random) and then change to a new i.i.d. 
value. We denote rji^k = l/^[7i,k] (^s the average change rate. 

This assumption is true for parameters such as user scheduling information, which are determined 
by BSs and stays constant per transmission time interval (TTI). On the other hand, parameters 



6 



such as channel fading may change continuously. However, the i.i.d. discrete time model on 
channel fading (named block fading model in other literature) is fairly accurate and widely used 



where M is the parameter. Note that E[7i,fc] = is unchanged under various values of M. This 
general distribution with different values of M can model different scenarios from deterministic 
process (i.e. Ti^k becomes deterministic as M — )■ oo) to Poisson process (i.e. TiM is exponentially 
distributed for M = 1). 

Because of the dynamics of these cooperation-dependent parameters, the end-user or its 
serving cell needs to detect the changes in their values. Several parameters like user scheduling 
information are determined by the BS and their values are thus known on real-time basis. For 
other parameters such as channel fading, the end-user can constantly measure their values through 
frequent pilot signals. Once the values of the parameters change, an overhead message will then 
be generated and sent. In this way, the overhead message will be updated every 7I,fc, a sufficient 
frequency without unnecessary burden on the overhead channel. 

After being generated by the serving base station BSi fc*, an overhead message is transmitted 
to BSi through inter-cell overhead channel. The overhead channel incurs delay denoted as 

fc. With this updating overhead message, BSj ^ can take the appropriate cooperation strategy 
accordingly. Note that each overhead message only has a lifetime of Ti,k, because the parameters 
change after Ti^k and a new overhead message is generated. 

We now consider overhead design in CoMP ZFBF, where coordinated neighboring BSs use 
zero-forcing precoders to null their interference |l2l, [fTTll . Therefore, the cooperation-dependent 
parameter in CoMP ZFBF is the channel direction information (CDI) 



where hj ^ is the A^^^ x 1 fading vector between BSj ^ and the end-user. In this paper, we assume 
uncorrected Rayleigh fading, i.e. each component of hj ^ is i.i.d. complex Gaussian CJ\f{0, 1). 
According to Assumption [U the fading hj ^ stays constant for a time 7^,^ and then changes 
independently, i.e. block fading [|32l . In this specific example, Ti^k is the channel coherence 
time. 



m, m. 

We assume a parameterized gamma distribution on 71,^ in Assumption [T] 




(2) 




(3) 
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In the beginning of each fading block, the end-user observes the new fading value hi ^. 
It searches through a codebook Cj ^ known by both itself and BSj ,t, which consists of 2^^ '= 
codewords. From the codebook Cj fc, the end-user will choose the codeword Cj ^ closest to current 
fading direction hj ^ such that |hjfcCjfc| is maximized. The index of this selected codeword Cj ^ 
is fed back to BSi fe* using i?j ^ bits. For BSi fc*, the overhead messages form an arrival process 
with inter-arrival time 7i,fc. The serving BS BSi fc* then transmits these overhead messages to 
BSj^fc through an inter-cell overhead channel. Based on the received overhead, BSj^fc chooses a 
zero-forcing precoding vector fj ^ such that |ft,A;Cj fcp = 0. See Fig. [T]for a conceptual plot of 
overhead messaging in CoMP ZFBF. 




Fig. 1. A conceptual plot of CoMP ZFBF in a heterogeneous cellular network. The end-user's serving BS coordinates with 
a pico BS. As the channel fading h between the pico BS and the end-user changes over time, the serving BS will perform 
frequent overhead messaging through backhaul, to inform the pico BS about the current fading values. 

C. The Impact of Overhead Delay 

Realistic overhead messaging has two major imperfections - delay and quantization inaccuracy. 
To make the discussion more concrete, we describe the impact of overhead delay in the context 
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of CoMP ZFBF, while the situations in other CoMP schemes are similar. The delay ^ of 
an overhead message is defined as the time between when that overhead is generated (i.e. the 
beginning of the respective fading block) and when it is received by BS/ fc- It is caused by 
unavoidable propagation time and the imperfections of the overhead channel such as congestion 
and hardware delays ll26l . We call this time window the overhead messaging phase. If the 
overhead delay Vi^k is smaller than the fading block length 71,^, we call the rest time Ti,k — T^i,k 
in that fading block the cooperation phase. Note that the cooperation phase may not exist in a 
fading block if the overhead delay Vi^k is larger than 71,^- See Fig. |2lfor an example. 
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Fig. 2. Overhead messaging pliases and cooperation phases of a coordinated BS in CoMP ZFBF. The overhead message phase 
starts from the beginning of each fading block and has a time length of overhead delay T). In a fading block, the coordinated 
BS will have cooperation phase only if the overhead delay V is smaller than the fading block length T. 

The interference from BS^ fc is different between these two phases. In the overhead messaging 
phase, the channel fading ^ has already changed but BSj ^ does not know its value yet. 
Therefore the zero-forcing precoder fj ^ is still determined from previously received overhead 
message. According to Assumption [TJ the current fading state hj ^ is independent of the previous 
fading block and thus the precoder ^ based on it. Therefore, statistically the interference 
|fj,fchjfcp is not reduced in the overhead message phase, which is the worst-case interference 
scenario [l23ll. [l24l. 

On the other hand, BSj ^ receives the new overhead message in the cooperation phase and 
adjusts its zero-forcing precoder fj jt accordingly. Because of Assumption [T} the fading value 
hj is assumed to keep unchanged during this block. Therefore the new overhead message 
remains accurate (minus quantization errors) in the cooperation phase and the selected precoder 
fj^fc minimizes the interference |fj_fchi,fcp for the entire phase ll23l . [|24l|. This is the best-case 
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interference scenario. 

It is seen that Assumption \T\ simplifies the impact of overhead delay into two opposite 
extreme cases: the interference to BSj ^ is either not reduced in the overhead messaging phase 
or maximally reduced in the cooperation phase. In practice, the channel fading hj ^ continuously 
changes with temporal correlation, instead of the i.i.d. discrete-time model in Assumption [U 
Therefore, the interference from BSj ^ should gradually change over time and is in fact bounded 
by the two extreme cases. However, it is hard to quantify mathematically. We use Assumption [T] 
because it is tractable for analysis, widely used in previous literature on MIMO systems [|23l . [|32l 
and allows a first-order analysis on the impact of overhead delay. Future work should consider 
more complicated models on channel fading, such as other discrete time models in [|25l . [|33l . 

m. 



D. The Impact of Overhead Quantization Error 

Another concern of realistic overhead messaging is the finite overhead quantization bits i?j fc. 
Larger Bj ^ usually translates into smaller quantization error and thus higher cooperation gains. 
However, the exact impact of Bj ^ depends on the specific CoMP scheme and the overhead 
codebook Cjjt. See [23] for an overview. We now discuss its impact in the context of CoMP 
ZFBF. 

The end-user's SIR 7 in CoMP ZFBF can be generally expressed ai3 

X,,feG U l-fcU^l.fc*} 

fc=i 

where the value of |fj fchj ^l is elaborated on in the following. 

1) The end-user's serving BS BSi fc* needs to null its interference to the L coordinated cells. 
Meanwhile it also wants to maximize the signal power |fi fc*hi fc*p to the end-user Thus 
its precoder fi fc* is chosen such that |fi fc*hi ~ X2Nk*-2L t^'^l- 

2) A coordinated BS BSj^^ G Bl cannot null its interference during the overhead messaging 
phase, because its zero-forcing precoder fj ^ is independent of hj^fc and we have |fj ^hj fcp ~ 



'in modem cellular networks, thermal noise is not an important consideration either in cell interior where the signal power 
is strong or in cell edge where interference is usually much larger. Interference is more dominant especially in HCNs because 
of additional overlaid cells with high spatial densities. We therefore neglect thermal noise and consider SIR in this paper. 
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X2- In the cooperation phase, BSj ^ receives the updated overhead indicating the codeword 
Cjfc = argmaXceCifc |ht,fcc| and chooses fj^ such that |fi,fcCjfcp = 0. However, because 
Cj is not the exact GDI hj ^, the value of |fj ^hj ^.p is still positive and dependent on the 
design of overhead codebook Cj fc. In this paper, we assume random vector quantization 
(RVQ) codebook d^k, which is commonly used in previous CoMP ZFBF [El, [IH, [l24l . 
Under this assumption, we then have Ifj^fchj^p ~ 2 ^fc-^ for BSj^ in the cooperation 
phase. 

3) A non-coordinated BS BSj ^ ^ Bl chooses the precoder independent of its interference to 

the end-user, and will have |fifchjfcp ~ xl- 
Based on the derivations above, the end-user SIR can be expressed as 

Yl ^ PkPi,kSi,k\^i,k\~'^'' 

where 5*1 fc* ~ x^{2Nk* — 2L), 5*^ ^ ~ ^'^d Pi,k is the interference cancellation factor for 

BSj^fc 

1 For BSi,fc i Bl 
Pi,k = \ I For BSjjfc G Bl during the overhead messaging phase (6) 

2 ^k-^ For BSj^fc £ Bl during the cooperation phase 

III. CoMP Throughput Evaluation With Imperfect Overhead Messaging 

In this section, we present a throughput evaluation framework for a class of downlink CoMP 
schemes, where inter-cell overhead messaging does not include user data. Common examples 
in this category include CoMP beamforming [fTTll . ||3TI and interference alignment [3J. The 
general framework considers the practical issues from overhead messaging but is also analytically 
tractable. 

Lemma 1: The long-term time fraction Tj ^ that a coordinated BS BSj ^ is in the cooperation 
phase is 

oo 

ri,k = p{Ti,k, oo) - r]i^k J {p{Ti,k, oo) - p{Ti,k, s)} ds, (7) 



where p{Ti,k,s) = P(r'i,fc < Ti,k,T^i,k < d) is overhead delay distribution. 
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Proof: See Appendix lAl ■ 
In our previous work [26], we provide general models on overhead messaging delay for both 
backhaul and over-the-air inter-cell overhead channels. We then derive the delay distribution 
p{-) as an explicit function of inter-cell overhead channel parameters (e.g. over-the-air overhead 
channel bandwidth). Note that our overhead delay models only affect the characterization of 
p(-), while the result in Lemma [H holds under various forms of delay distribution p{-). 

In previous literature, no overhead delay is considered, i.e. delay Vi k = and thus p{Ti^k-, s) = 
1 for any s and 7i,fc. Under this condition, ^ in (|7]) becomes 1 as well, which means BSj^fc 
is always in the cooperation phase with interference maximally reduced. This is obviously an 
over-optimistic prediction on the interference. 

Theorem 1: For CoMP schemes without user data sharing, the end-user's long-term throughput 

is 

n=J2 pbR{ib), (8) 

where the summation is over all possible subset B C Bl, 'Jb is the end-user SIR when BSs G B 
are in the cooperation phase and BSs e Bl \ B axe in the overhead messaging phase, R{-) is 
the SIR-rate mapping function, and 

I n ^^fc n (i-^^-fe) I ' (9) 

Proof: Due to the overhead messaging delay, each coordinated BS BSj^^ G Bl now has two 
states: 1) the overhead messaging phase (with probability 1 — fc) when the overhead message 
at BSj fe is already outdated and the updated overhead has not been received yet; and 2) the 
cooperation phase (with probability Ti^k) when BSj ^ receives the updated overhead message. 
BSj^fc has different cooperation performance between these two states, for example, as shown in 
^ for CoMP ZFBF. Therefore, we use a subset B C Bl to denote the scenario that only BSs 
G B are in the cooperation phase. The probability of this scenario is 

PB = I n n (1 - ^^^) I • (10) 

\BS,,fceB BS,„keBL\B J 

Each subset B C Bl corresponds to a possible scenario regarding which coordinated BSs are in 
the cooperation phase. The end-user's long-term is the average rate over all possible scenarios. 
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As shown in ([8]), the CoMP throughput evaluation framework explicitly quantifies the impact 
of overhead delay through the time fraction Tj ^ and probability distribution p^. It also includes 
the impact of finite overhead bit size because the end-user SIR 73 is affected by overhead 
quantization error, for example, as shown in Q for CoMP ZFBF. In this way, the framework 
considers the imperfections in the overhead messaging. Combined with the SIR characterizations, 
it can be used to quantify the throughput and coverage of different CoMP schemes. 

As we mentioned before, previous work does not fully consider practical issues in inter-cell 
overhead messaging. In particular, they ignore the overhead delay, which is in fact non-trivial in 
most network environments |[9l, ifTTI . lfT3l . To show the importance of modeling and analysing 
overhead delay in CoMP evaluation, we compare the result in Theorem \T\ with previous work. 
Without overhead delay, the coordinated BSs always have the updated overhead and stay in the 
cooperation phase. Obviously, this is a idealized special case of Theorem [U 

Corollary 1: Under the assumption of delay-free overhead messaging, the long-term CoMP 
throughput is 

7^ = i^(7^3J, (11) 

where 75^ is defined in Theorem \T\ 

Proof: When overhead messaging has no delay, all coordinated BSs will always be in the 
cooperation phase, i.e. Tjjt = 1 for each BSj^fe G Bl- Therefore we have ps = 1 if B = Bl and 
Pb = otherwise. The end-user's rate is then TZ = Ri^Q^). ■ 
As shown in Corollary [H assuming no overhead delay significantly simplifies the analysis 
of CoMP schemes. However, as we elaborated before and will be shown in the numerical 
simulations, this assumption is far from the reality and causes the wide gaps between analytical 
predictions and realistic implementations. 

IV. CoMP Zero-Forcing Beamforming Throughput Analysis 

In this section, we derives the distribution of the end-user's SIR 7^ for CoMP ZFBF, which will 
be used with the evaluation framework in Theorem [T] to quantify the CoMP ZFBF throughput. 
We first derive the SIR CDF in 1-tier cellular networks in Theorem |2l and then extend it to the 
general HCN scenario in Theorem [3l 

For now, we simplify the notation for 1-tier cellular networks by dropping the tier index k. 
Specifically in 1-tier networks, the BSs have the same transmit power P, number of antennas 



13 



and path loss exponent a. Their locations {Xj, i G N} form a PPP $ with intensity A. Therefore 
the serving BS is simply BSi - that is, the nearest BS to the end user - and the L coordinated 
BSs are = {BS2, . . . , BS^+i} - that is, the second to (L + 1) nearest BSs to the end-user. 
Similarly, we now use notation p,; and Bi as simplified versions of p,; ^ and Bi ^ in 
In 1-tier networks, the end-user SIR 7^ in ([8]) can be simplified as 



PSi 




—a 


Si 


Xl 


—a 


E PP^S^\X, 

Xie^\{Xi} 


X,G*\{Xi} 


—a 



where pi = 2^-1 for BSj G B and pi = 1 otherwise. 

Theorem 2: In 1-tier cellular networks, the CDF of the end-user SIR 7^ in 



(12) 



< 



/3 



N-L- 



i=2 



r(i+f)(^-i)! 



r ^^ 



> l-exp|-r(l + 2/a) 
where / = is the cardinality of the set B, p 



{N-L)T{l-a/2) 



/3[3-°'pmin + (2«+3)-'»(l-pmin)] 



is bounded as 



(13) 



mm = min Ipi}, and 



1 



pi 



BSi ^ S 

2-1^ BSi G i3 
Proof: We denote the normalized interference (normalized by P) as 

^,e*\{Xi} BSigB BS,e<i>\{BSiU'B} 

The last equality comes from the definition of B. 

Upper Bound: For a PPP $ = {Xi,X2, . . .} inlR^ with an arbitrary density A, {|Xip, jXg 

form a one dimensional PPP with intensity ttA. We thus have 



(14) 



(15) 



E 



|Xi|" 

IX; h 



E 



1^ 

IX; I 



2 \ a/2' 



r(i + f) (z + 1)! 

r(z + f) 



(16) 



where the last equality holds from Appendix |Bl The upper bound on the SIR CDF is 



P 7 



5i|Xi|-° 



</3 



P 



> 



1\ (") 



5i|Xi|-" - /3 



< /3E 



^1 



-^el-^il 



/3E 



1 



5^p,E[5,]E 



i=2 



(17) 

1 



N-L-l 



where (a) follows from Markov's inequality. As elaborated below ([5]), we have E 
(because ~ X2N-2l) ^i^^ W\S.-] = 1 (because Si ~ X2) i > 2. Therefore Combing ([16 
and (fT71) gives the CDF upper bound in Theorem [2l 
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Lower Bound: Let /(m) = ^Xie<i>\{Xi x^} other words, /(m) is the sum 

interference experienced by the end-user, if the nearest m BSs (including the serving BS) are 
removed or turned off. Note that /(q) means that the end-user is not associated with any BS, and 
experiences interference from all the BSs in the PPP $. Denote / = \B\ as the cardinality of B 
and pmin = ininBSieelPi}- The actual interference can be expressed as a function of /(„) 

PminSi\Xi\ " + >S'i|Xj| " 

BS.eB x,e<s>\{Xi(jB} 

BS,G{BS2,...,BS,+i} X,g*\{Xi,...,X,+i} 

= Pmin 'S'j|Xj| " + (1 — Pmin) 'S'jIXjl " 

X,&'i>\{Xi} X,G$\{Xi,...,X,+i} 

= Pminl(l) + (1 — Piain)I(l+l) 

h + (2/ + 3)-"(l - p^in)] /(O) (18) 

where ^ means stochastic dominance. In the right hand side of (b), we assume the strongest / 
interfering BSs are cancelling their interference, instead of the / BSs in the set of B. Therefore 
it is a lower bound on the interference Jg. (c) holds from Appendix O deriving the lower bound 
of an arbitrary /(„). 

Based on the lower bound on interference J, the SIR CDF can be lower bounded as follows 

V Ib - - V 3-"Pn.i. + 2/ + 3 -'^y - ^ J 



(19) 



^h^^^ ^ = m--f^^Stl^)\i-P^..)V ^(^) = + 1)^-^, and (d) holds from Theorem 1 

in [|35l . The CDF lower bound is then derived as 

H < /3) > 1 - E [e-^(^)] 

> 1 - exp |-7rAr(l + 2/a)(E[Z])^^| 

(iV-L)r(l -a/2) T"'/"' 



= l-exp|-r(l + 2/a) 
where (e) holds from Jensen's inequality. 



/3[3-Vmin+(2/ + 3)-"(l-pmin)]. 



(20) 
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Based on the results in Theorem [2] for the 1-tier networks, we now characterize the end-user's 
SIR 7e in a general i^-tier HCN. 

Theorem 3: In a general A'-tier heterogeneous cellular network, the CDF of the end-user SIR 
7/3 in dH]) is bounded as 



< 



JV,,*-L-1 



oo r(i+^)(^-i)! ^ p^, - (A,..)-^r(i+T)(--i)! 

^'«,A:* r/'i^^'l "T Pi.* ^ ^^-'^ 

fc=l j=l 



1=2 



> 1 — exp 



ttA 



amax p / _|_ 2 



(jv,*-L)r(i-^) 



/3[3-«maxp^i„ + (2i+3)-<^max(l_p„^i„) 



(21) 



K 



where / = X = ^ {Pk/Pk*)"'^, "max 
fe=i 



max{ai, . . . , ok}, Pmin = min {pi,k}, and 



1 



BSi,k i B 



(22) 



2 BS,,fe Gi3 

Proof: In the i^-tier HCN, the normalized interference Jg can be written as 

^ K p 

where pj^fc = 2'^'=-^ for BSjjt G i^z, and pi^k = 1 otherwise. 

Upper Bound. Similar to the proof for the 1-tier case, the upper bound of SIR CDF is 

K oo ^ 



(23) 



P 7 



Si,k* 1^1,. 



</3 < 



/3 



N-L-1 



1=2 



1^1, 



k=l i=l 
k^k* 



Following the same steps in Appendix |B] we have 



E 



1^1, 



IX- 



i,k* 



r(i + ^) (2-1)! 
r(. + ^) 



For k 7^ k*, Xi^k and Xi fc* belong to two independent PPPs. Therefore we have 



E 



1^1. 



E [|Xi,fc*r^*] _ (Afc*7r)-Tr (1 + 1)! 



EflX.|°^l 



(Afc7r)^r(^ + f) 



|Xj^fc 
(24) 

(25) 
(26) 



Therefore the upper bound in Theorem [3] is proven. 

Lower Bound. The key idea of this proof is converting the interference from K tiers to that 
of a single tier. Then we apply the lower bound from Theorem [21 Comparing (|23|) with (fT5l) . 
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it is seen that the i^-tier case is different from the 1-tier case in two important aspects: 1) BSs 
from different tiers have different powers; and 2) BSs from different tiers have different path 
loss exponents. 

We first present the way of eliminating the power differences. The interference from the /c*'* 
tier is 



Ik — PkPi,kSi^k\Xi^k\ — PkPi,kSi^k 



(Pk/Pk*)"^ 



0!k 



^ ^ Pk Pi,kSi^k\Yi^k\ 

(27) 



where Yi^ = — j_ . The conservation property in ||35l . Il36l states that {Yik,i G N} form 

a new PPP with intensity = (Pk/Pk*)"'' ■ Therefore the interference can be viewed as 
generated from the K new tiers {l>i, . . . , <l>x} with the same transmitting power Pk*. Therefore, 
the normalized interference in (|23|) can be rewritten as 

K 

^^ = 11 H P^,kS^,k\YiA"'^ (28) 

where = means equivalence in distribution. 

We then set all the path loss exponents to a common value «max = max(ai, . . . , uk), which 



is the best case since the interference attenuates faster. The normalized interference in (1231) can 
be lower bounded as 

K 

Ib>Y, H P^,fc5,,fc|r,,fc|-"— = 4^ (29) 

K ^ 

Denote l> = |J <lfc. Because {$1, . . . , l>i^} are independent PPPs, $ is also a PPP with 

fc=i 

intensity A = Yl, ^k- The interference lower bound can be viewed as generated from a single 

k=l 

tier where 1) the BS locations forms a PPP $ with intensity A; 2) BSs have the same transmitting 
power Pfc*; and 3) BSs have the same path loss exponent amax- Therefore, can be rewritten 



as 



— / ^ Pi,k^i,k\J^ i,k\ 

^.,fee*\{^i,fc*} 

= Pi,kSi,k\Yi,k\ + Si^k\Yi^k\ (30) 
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It is obvious now that /g is in the same form as Jg in (fT5l) . 

Applying (fT9l) and (|20|) . we have the lower bound on the end-user SIR 



P(7 = -^'-'^-^'-^^ </3) >p( -^^-'r^.^k*^ <^ 



> 1 — exp 



-vrAr ( 1 + ) (E[Z])-^ 



(31) 



where Z = ;3[3-^n.axp!;,'+(f/;'3)LI. (i,^,^^,^] , l = \B\, and p^i^ = minBs,,eB{A,fc}- From Lemma 

1 and 3 in Wl, we have E[|Xi,fe*|-"''*] = (7rA)Tr(l - ak*/2), which gives 

(iVfc* -L)(7rA)Tr(l- 
]g 2^ = ^ ^ V 2_^! ('32') 

Therefore the lower bound in Theorem [3] follows. 

■ 

The bounds in Theorem [2] and [3] are insightful as they show clear dependence on important 
parameters such as the number of BS antennas N, the path loss exponent a and the number 
of coordinated neighbouring BSs L. On the other hand, the bounds in 1-tier network case are 
independent of the BS spatial density A. This is often called scale-invariance which is a known 
property of interference-limited cellular networks EOl . [|27l . Il37l . These bounds on the end-user 
SIR 70 can be used with the throughput evaluation framework in Theorem [T] to derive the bounds 
on the end-user throughput. 

V. Numerical Results and Discussion 

In this section, we simulate the CoMP ZFBF performance under both realistic and idealized 
overhead messaging. We consider a 3-tier heterogeneous cellular network comprising macro 
(tier 1), pico (tier 2) and femto (tier 3) BSs. We simulate the scenario that the serving BS is a 
macrocell BS, i.e. k* = 1. Notation and system parameters are given in Table HI 

Two main performance metrics for CoMP schemes (including CoMP ZFBF) are their improve- 
ments on network coverage and capacity. We simulate both performance metrics, by considering 
the different types of SIR-rate mapping functions R{-) in Theorem [B 

1) (Coverage under CoMP ZFBF) If the end-user only requires a fixed target rate T^target 
(e.g. a voice user), its SIR-rate mapping function R{-) is 

(33) 

I T < Hnrgtt 
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TABLE I 

NOTATION & Simulation Summary 



notation 


Description 


Simulation Value 


K 


The number of tiers in the HCN 


3 


k 


The tier index 


k=l (macro), 2 (pico), 3 (femto) 


k* 


The tier index of the serving BS 


1 (macro) 


Pk 


Transmitting powers of fc*''-th tier 


Pi = mW, Pa = 2W, Pa = 0.2W 


Xk 


Spatial density of fc*'' tier 


Ai = 10-V"^^ A2 = 10-V^n^ A3 = IQ-^/m^ 




Path loss exponent of fc"* tier 


Qfi — 4.0, 02 = 3.5, 0:3 — 3.0 


Nk 


The number of BS antennas in the fc-th tier 


A^i = 8, iV2 = 4, A^a = 2 




Overhead delay 


Not fixed 




Time fraction of BSi,fc in cooperation phase 


Not fixed 


Ti,k 


Channel fading block length of BSi,fc 


Fixed length of 80 ms 


Bi,k 


Overhead message quantization bits for BSi^fc 


Not fixed 


L 


The number of coordinated cells 


Not fixed 


Bl 


The set of coordinated cells 


Not fixed 


7target 


The target SIR used in ([33} 


3 dB 


G 


The Shannon gap used in ([34} 


3 dB 



The throughput quantified in Theorem [U is then simply the user's target rate times its 
probability of being in coverage (i.e. with SIR 7 larger than the target SIR 7target)- 
Therefore, we can quantify the coverage improvement from CoMP ZFBF by normalizing 
the derived throughput by T^target- We use 7target = 3 dB in the simulations. 
2) (Throughput under CoMP beamforming) If the end-user is a data-greedy user, its SIR-rate 
mapping function R{-) is 

^(7) = log2 (1 + ^) bps/Hz, (34) 

where G is the Shannon gap, which is 3 dB in our simulations. Quantifying the rate of 
users of this kind will show the throughput improvement from CoMP ZFBF. 

A. The Configurations of Overhead Channel 

Regarding inter-cell overhead channels, we consider the scenario that the overhead messages 
are shared through the BSs' backhaul. In our previous work [|26l . the backhaul connection 
between two coordinated BSs is modelled as a tandem queue network consisting of several 
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servers (e.g. switches, routers and gateways), each of which has exponential processing time. 
See [26] for more details on this model. The limited processing rate from the backhaul servers 
inevitably introduces overhead delay, whose distribution piTi^k, d) is derived in Theorem 1 in 
||26l and used in the simulations. 



0.82 




20 40 60 80 100 120 140 160 

The average delay (ms) 



Fig. 3. Downlink CoMP ZFBF coverage probability vs. the average overhead channel delay. The overhead bit size is Bi^k = 
3(A'^fc — 1), which gives pi^k ~ 12.5% (i.e. a coordinated BSi,fc can cancel 87.5% of its interference once receiving the updated 
overhead). The number of coordinated cells is L = 1. 



O CoMP ZFBF with L=1 under realistic overhead messaging 
—A— CoMP ZFBF with L=1 under delay-free overhead messaging 
3.9 -B- No CoMP ZFBF 




20 40 60 80 100 120 140 160 

The average delay (ms) 



Fig. 4. Downlink CoMP ZFBF throughput vs. the average overhead channel delay. The values of L and Bi.k are the same as 
Fig. [3] (i.e. B,,k = 3(iVfc - 1) and L = 1). 
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Fig. [3] and |4] show the impact of backhaul overhead delay on the CoMP ZFBF coverage 
and throughput. It is seen that CoMP ZFBF coverage and throughput fall almost linearly as 
the average delay grows from zero (i.e. delay-free overhead channel as assumed in previous 
literature). In fact, when the overhead channel delay is larger than 60% of the fading coherence 
time 7i,fc5 CoMP ZFBF does not bring any coverage or throughput gain. This observation diverges 
significantly from previous optimistic performance predictions which ignore the overhead delay. 
It also provides a rule of thumb for the overhead channel configurations for CoMP ZFBF. 

B. Choosing Coordinated Cells 



0.85 




1 2 3 4 5 

The number of coordinated BS 



Fig. 5. CoMP ZFBF coverage probability vs. the number of coordinated cells L. We use Bi,k — 3{Nk — l) to give pi^t ~ 12.5%, 
i.e. a coordinated BS BSi.k can cancel 87.5% of its interference once receiving the overhead. For the overhead delay Di^k between 
the serving cell and BS^^fe, we adjust the servers' processing rates in their backhaul path, to make sure that the average overhead 
delay E[I?i,fc] = 20 ms fU]. 

A fundamental design question for CoMP beamforming is how many and which neighbouring 
cells should be coordinated. Coordinating more cells may translates into less interference from 
other cells, but also weaker signal power for the end-user and heavier overhead messaging 
burden. Fig. \5\ and |6] show that the number of coordinated cells should be kept fairly small in 
HCNs, even under ideal overhead model (i.e. no overhead delay and infinite quantization bits) 
and limited feedback overhead model (i.e. no overhead delay but finite quantization bits). This 
observation diverges significantly from previous work in conventional macrocell networks and 
implies that dominant interference comes from only a few neighboring cells in HCNs. Further, 
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1 2 3 4 5 

The number of coordinated BS 



Fig. 6. CoMP ZFBF throughput (bps/Hz) vs. the number of coordinated cells L. The configurations on Bi^k and the backhaul 
channel are the same as Fig. ID (i.e. B,,fc = 3{Nk - 1) and E[D,,fc] = 20 ms). 

by considering both overhead delay and rate constraints, our work indicates that coordinating 
with only one other cell (i.e L = 1) is actually optimal for a serving base station BSj fc* with 
eight antennas. 



g 

'm 14 




jC 20 40 60 80 100 120 140 160 

The average delay (ms) 



Fig. 7. CoMP ZFBF coverage and throughput losses (in percentage) under intra-tier coordination vs. the average overhead 
channel delay. The values of L and Bi,k are the same as Fig. |3] (i.e. Bi^k ~ 3{Nk — 1) and L = 1). 



When the serving base station BSi^a,* coordinates with L other cells, we assume the optimal 
scenario that the L strongest interfering cells (possibly coming from different tiers) are able 
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to coordinate with the serving cell. In practice, cross-tier coordination may be restricted in 
HCNs and only intra-tier coordination is allowed. For example, the end-user installed femtocells 
are controlled by their owners, and may not be allowed to or capable of coordinating with 
macrocells [[T6l . Coordinating only with cells in the same tier (termed intra-tier coordination) is 
of course sub-optimal and results in coverage and throughput losses, which are quantified in Fig. 
|7](w.r.t. different overhead delay profile) and[8](w.r.t. different values of L) for CoMP ZFBF The 
performance loss is most significant (e.g. as high as 20%) in the practically important scenario 
of CoMP ZFBF - that is, relatively large number of coordinated cells and/or small overhead 
channel delay (e.g. if the overhead channel is optimized as in [[TOl '). 
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Fig. 8. CoMP ZFBF coverage and throughput losses (in percentage) under intra-tier coordination vs. the number of coordinated 
cells L. The configurations on Bi^k and the backhaul channel are the same as Fig. [5] (i.e. Bi^t = 3(A'^fe — 1) and E[I'i.fc] = 20 
ms). 



VT Conclusion 

This paper presents a novel approach to evaluate downlink CoMP schemes in the new network 
paradigm of HCNs, by developing a new framework to quantify the impact of overhead delay 
and using PPP model in end-user SINR characterization. This proposed approach can be used 
for a class of CoMP schemes, and is applied to CoMP ZFBF in this paper as an example. 
We show that CoMP ZFBF performance heavily depends on the overhead delay and its design 
should be fairly conservative (e.g. coordinating with only one or two other cells). These results 
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align with the findings from several industrial implementations ll9l- lfTT1l . and provide insights 
diverging significantly from previous work without overhead delay modeling. 

In this paper, one fundamental assumption is the i.i.d. discrete time model of the cooperation- 
dependent parameters. Developing more complicated and accurate models is an important related 
topic for future works. Besides, the evaluation framework assumes the number of coordinated 
cells L as a pre-determined parameter. The adaptive selection on L as a function of instantaneous 
other-cell interference can potentially bring more CoMP gains, and should be considered in future 
CoMP evaluation in HCNs. Considering multiple user antennas in CoMP is also a related topic 
for future research. 



A. Proof of Lemma \T\ 

The cooperation phase only occurs in fading blocks for which the overhead messaging delay 
V is smaller than T. Here we omit the subscripts of delay V and block length T to keep the 
proof general. The percentage of these fading blocks is 



where the last equality holds by definition. In these fading blocks, the overhead messaging phase 
will have a time length E[I?|I? < 7] 



Appendix 



P(p < T) = ¥{V <r,V < oo) = p{T, oo) 



(35) 



oo 




(36) 







Thus the average duration of the cooperation phase is IE[7~] —E,[D\V < T]. 
In sum, the long-term time fraction of the cooperation phase is 




oo 



p(r, oo) 



p{T, oo) 

nr] 



{1 - F{V < s\V < r)}ds 







oo 




(37) 







By definition, p(r, oo) = F{V <T,V <oo) = F{V < T). Therefore (a) follows. 
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B. Auxiliary Result for the CDF Upper Bound in Theorem.^ 

Consider a one-dimensional PPP $ = {Fi, F2, • • •} with intensity A, we have 

i 

where Aj ~ exp(A). For an arbitrary positive number u, the following equality holds 



(38) 



E 









= E 


[( 









Ai 



E 



Ai + A2 + . . . A, 

r(i + z/) u-i)\ 



E <^ E 



Ai 



Ai + A2 + . . . Ai 

r(i + z/) {t-iy. 

T(i + u) 



Ai + . . . Ai = X 



C. Auxiliary Result for the CDF Lower Bound in Theorem |2] 
For an arbitrary m, /(m) can be expressed as 

X,€<S>\{Xi,...,X^} 
= Si (|Xj| — \Xjn\ + 

X,€<S>\{Xl,...,Xm} 



= E >i 

x,e*\{Xi,...,x™} 

X,G'I>\{Xi,...,X™} 



1 + 





Xm 













Xm 1 




1 Y 


1 







Si{\Xi\ — \Xr 

—a 

S^{\X,\ - 



x„ 





Xm\ 




1 Y 


1 




1 Y 1 



5^ S,{\Xi\ - 

X,G<I'\{Xi,...,X„} 



Since the above inequality holds for any realization of \X„i\ and we have 

E[|X„|] 



> 1 + 



E[|X,„+i|]-E[|X,„|] 



s,{\x,\ - \x^\y 

x,e$\{Xi,...,x„} 



(a) 



:i + 2m)-° S,(\X,\-\X^\)-\ 

X,G*\{Xi,...,X„} 



where (a) follows because we have E[|Xm|] = (Att 

A 



0.5 r(m+0.5) 
' (m-1)! ■ 



(39) 



(40) 



(41) 



Now we define another point process = {Yj E : for any, Yj = Xi — = 
i — m}, i.e. $ is formed by moving every point of {X.m+iiXm+2^ ■ ■ ■} in the PPP $ toward 
the origin by distance \Xm\- Note that is also a spatial Poisson Point Process, but non- 
homogeneous with a larger density than the original $ (because when moving points toward the 
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origin, we actually compress the space in R^). Therefore the sum interference from is larger 
than that from the original i.e. larger than /(o). 

oo 

hm) >■ (1 + 2m)-" - 

i=m+l 

oo 

= (l + 2m)-" Yl Sj+m\Yj\-^ 

b (l + 2m)-"/(o). (42) 
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